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Abstract 

In order to mimic the phase changes in the primordial Big Bang, several cosmological solid-state 
experiments have been conceived, during the last decade, to investigate the spontaneous symmetry 
breaking in superconductors and superfluids cooled through their transition temperature. In one of 
such experiments the number of magnetic flux quanta spontaneously trapped in a superconducting 
loop was measured by means of a long Josephson tunnel junction built on top of the loop itself. 
We have analyzed this system and found a number of interesting features not occurring in the 
conventional case with simply connected electrodes. In particular, the fluxoid quantization results 
in a frustration of the Josephson phase, which, in turn, reduces the junction critical current. 
Further, the possible stable states of the system are obtained by a self-consistent application of 
the principle of minimum energy. The theoretical findings are supported by measurements on a 
number of samples having different geometrical configuration. The experiments demonstrate that 
a very large signal-to-noise ratio can be achieved in the flux quanta detection. 

PACS numbers: 74.50.+r, 85.25.Cp, 98.80.Bp 
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a) Sketch of an in-line Josephson tunnel junction with a single doubly con- 
nected ring-shaped base electrode. Lioop = Li + L2 + Lf,. b) Sketch of an 
in-line Josephson tunnel junction with two doubly connected electrodes. The 
base electrode is in black, the top electrode is in gray, and the tunneling in- 
sulating layer is white. The dotted withe arrows indicate the direction of the 
circulating currents [5] 

(a) Top view (not to scale) of an in-line long Josephson tunnel junction in 
the most general bias configuration, (b) Cross section (not to scale) of the 
junction area showing the currents in the electrodes. The base electrode is in 
black, the top electrode in gray and the tunneling insulating layer is white in 
(a) while it has a wavy pattern in (b) [9] 

Elementary cell of the equivalent lumped circuit for the static Josephson 
transmission line [TT] 

(a) Phase profile 0(x) of a long {i = 10) in-line symmetric {a = 1) Josephson 



tunnel junction as in Eq.(17) for hr = —hi = 2 that is ^i^r = 0; (b) its first 
derivative (j)x, (c) its second derivative (pxx, and (d) its normalized energy 
density Tisc 

Numerically computed zero- field normalized junction energy Hsg vs. l"^ of a 
symmetric {2aAb = 1) in-line Josephson tunnel junctions with for different 



normalized lengths i. For £ — )■ 00 the analytical expression in Eq.(23) was 



plotted. L = I/Iq is the normalized bias current [T7] 

6 Magnetic diffraction patterns ic{he) for a very long in-line Josephson junction 
with different values of the symmetry parameter < a < 1 and for three 
Afc values, namely, 0.25, 0.5 and 0.75. The critical current is normalized to 
JcWXj and the magnetic field to JcXj UHl 

7 Geometrical layouts (not to scale) of our samples: (a) rectangular and (b) 
circular loop. The semi-angular length of the junction is 5 ~ 1 |26] 



8 Experimental magnetic diffraction patterns of a long Josephson tunnel junc- 
tion built on top of a rectangular loop (sample JJB) under different bias 
configurations and applied magnetic field orientation: (a) in-plane and (b) 
transverse magnetic field. With reference to Fig. [7]^b), the voltage is mea- 
sured across the terminals 1 and 2, while the bias current is applied through 
the contacts: 1 and 2 for curves I, 1 and 3 for curves II and 1 and 5 for curves 

III. Jo = 460 /iA, Hy,c = 490 A/m, and Hz,c = 80 A/m at T = 4.2 K |27] 

9 Variation of the critical current as a function of the magnetic flux with which 
the loop was field-cooled through the superconducting transition temperature. 

Hy = Hz = and T = 4:.2K M 



I. INTRODUCTION 

Long Josephson tunnel junctions (LJTJs) were traditionally used to investigate the 
physics of non-linear phenomeneP^. In the last decade they have been employed to shed 
light on other fundamental concepts in physics such as the symmetry principles and how 
they are brokerP^I. A recent experimeniP has demonstrated spontaneous symmetry breaking 
during the superconducting phase transition of a metal ring and both fluxoids or antiflux- 
oids can be trapped in the ring while it is cooled rapidly through the superconducting 
critical temperature. The basic phenomenon of quantization of magnetic flux in a multi- 
ply connected superconductor was suggested long time ago as one among several possible 
condensed matter cosmological experiment^ suitable to check the validity of the causality 
principle in the early Univers^. In the experiment of Ref.^ the magnetic flux quanta are 
spontaneously trapped in the ring during its cooling through the transition temperature. 
Much later at lower temperature when superconductivity is fully established, the number 
0, ±1, ±2 ... of flux quanta is registered as a function of the quench rate. This can be done 
in a variety of ways, one of which is the detection of the induced persistent currents by the 
magnetic fleld modulation of the critical current of a planar LJTJ built on top of the ring. 
In these experiments, the quench rate can be varied over four decades. This allows for an 
accurate check of the theoretical predictions of the involved second-order phase transitions. 
This is of interest within cosmology and of major importance for the physical understanding 
of many order-disorder processes. However, the working principles of that experiment had 
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FIG. 1. a) Sketch of an in-line Josephson tunnel junction with a single doubly connected ring- 
shaped base electrode. Ljoop = Li + L2 + -^6- b) Sketch of an in-line Josephson tunnel junction 
with two doubly connected electrodes. The base electrode is in black, the top electrode is in gray, 
and the tunneling insulating layer is white. The dotted withe arrows indicate the direction of the 
circulating currents. 

not yet been reported. The general task of this work is to study the static properties of a 
planar LJTJ for which at least one of the superconducting electrodes is multiply-connected, 
i.e., not every closed path can be transformed into a point. In the simplest case, one of the 
superconducting thin-film stripes forming the LJTJ is shaped as a doubly-connected loop. 
This configuration is illustrated in FigQa) in which the ring-shaped base electrode is in 
black, while the top electrode is in gray and the junction area is in white. The geometry of 
the loop is not critical to our discussion; however, a ring-shaped bottom electrode simplifies 
the analysis. 

For the sake of generality, in our analysis, we will include an external flux $e linked to the 
loop by some externally applied field and the presence of an integer number n of flux quanta 
trapped in the loop; altogether they induce a current Idr = {n^o ~ ^e)/Lioop circulating 
clockwise in the loop and inversely proportional to its inductance, Lioop', in turns, the circu- 
lating current produces at the loop surface a radial magnetic field, Hrad oc Idr, that adds 
to any external field, Happ, applied in the loop plane. With no loss of generality, we will 
assume that the width, Wt, of the top film does not exceed the width, Wh, of the bottom 
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one, Wt < Wb and, to simplify the analysis, we will also assume that both widths are much 
smaller than the mean radius, R, of the ring; in this narrow ring approximation, the current 
distribution in the ring and the surface radial field become radially independentP. A dc cur- 
rent I is injected into the loop at an arbitrary point O along the ring and is inductively split 
in the two loop arms before going through the LJTJ; let a (1-a) be the fraction of the bias 
current I diverted in the left (right) side of the loop. In principle, a values outside the [0, 1] 
range are possible if the current / would include also a persistent current Idr circulating in 
the loop: however, since the two currents are independent, we will treat them separately. 
Independently of the a value, the bias current I is extracted at one end of the junction via 
the top electrode. With the current entering and exiting at the junction extremities we have 
the well-known case of the so-called in-line configuration treated in the pioneering works on 
LJTJs soon after the discovery of the Josephson effeclP'^. 

Throughout the paper we will limit our interest to LJTJs in the zero-voltage time- 
independent state; this can be achieved as far as the applied current I is smaller than 
the junction critical current Ic- To further simplify the analysis, we assume that the Joseph- 
son current density Jc is uniform over the barrier area and that the junction width W 
is smaller than the Josephson penetration depth Aj = a/ ^0/211 j^Qd^Jc setting the length 
unit of the physical processes occurring in the Josephson junction (here $0 is the magnetic 
flux quantum, /xq the vacuum permeability and de the junction magnetic thickness). The 
gauge-invariant phase difference (p of the order parameters of the superconductors on each 
side of the tunnel barrier obeys the Josephson equation^^^ 

Jz(X) = J,sin0(X), (1) 

and 

KV(j){X) = H X n, (2) 

in which — L/2 <X<L/2isa curvilinear coordinate and L is the long dimension of the 
junction. The net current crossing the tunnel barrier is I = W j^^i^Jz{X)dX . The last 
equation states that the phase gradient is everywhere proportional to the local magnetic field 
H and parallel to the barrier plane. Therefore, in the case of a curvilinear one-dimensional 
junction, a uniform external field applied in the junction plane has to be replaced by its radial 
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component!^, k = $o/27r(ie/io = Jc^'j has the dimension of a current {k ^ 2.5 mA when 
de ~ 100 nm, which is typical of all-Niobium Josephson junctions) and fi is the versor normal 
to the insulating barrier separating the two superconducting electrodes. It is well knowrP^ 
that combining Eqs.Q and ^ with the static Maxwell's equations, a static sine-Gordon 
equation is obtained that describe the behavior of a one-dimensional in-line LJTJ: 

A^,0=sin</,(X). (3) 

Equation([3]) was first introduced in the analysis of asymmetric in-line LJTJs by Ferrel and 
Prange^ in 1963; few years later, Owen and Scalapind^^ reported an extensive study of its 
analytical solutions for symmetric in-line Josephson junctions (provided that L > 7cXj/2). 
Ampere's law applied along the barrier perimeter requires that the magnetic fields at the 
two ends of the junctions differ by the amount of the enclosed current: / = W[Hy{Ij/2) — 
ify(— L/2)]. We remark that Eqs.Q, ^ and ^ automatically satisfy the Ampere's law. 
As it is usually done in the modeling of Josephson interferometers, it is useful to divide 
the loop inductance Lioop in three inductive paths characterized by positive coefficients Li, 
L2 and Lh having units of inductance such that Lioop = Li + L2 + so that any current 
Icir circulating around the loop sees them in series; furthermore, a = 1 (a = 0) in the 
limit of L2 » Li + Lb {L2 « Li + Lf,). More specifically, Li^2 is the angular fraction 
'i9i,2/27r of the inductance for an isolated superconducting narrow (Wb « R) ring* Lrmg = 
/xo-R(ln 8R/Wb—2+ln 4) and Lb, that should not be mistaken as the inductance of the LJTJ, is 
given by the junction physical length L times the inductance per unit length of the junction 
bottom strip-line Cb, related to the magnetic energy stored within a London penetration 
distance of its surface. If the bottom and top superconducting films have, respectively, 
thickness db and dt and bulk London penetration depths Xib and Xn, then, in presence of a 
quasi-static magnetic field, their effective penetration depthP^is Xb^t = XLb,ttainh{db,t/'^XLb,t) 
which reduces to XLb,t in the case of thick superconducting films {db,t > 5XLb,t)- The magnetic 
penetration de of a tunnel barrier with negligible height tox << db^t df, ^ Xb + Xf. 
Insofar as the width Wb is much largei'^^'^ than the strip-line magnetic thickness d^, then 
Cb — ^J'oXb/Wb^, this expression also takes into account the kinetic inductance due to the 
motion of the superelectrons. In the wide-strip approximation, most of the magnetic energy 
is confined in the region between the plates and the fringing field can be ignored; as the strip 
width becomes narrower, the fringe field effects become more important and may dominate 

7 



if Wb and de are comparabl^^. It is worth pointing out that, since, in all practical cases, the 
width of the loop is much larger than the London penetration depth, then is considerably 
smaller than the inductance per unit length along the ring Lnng/'^T^R] this is due to the 
presence of a counter electrode acting as a superconducting ground pland^^'^. Similarly, we 
introduce Ct — ^oXt/Wt, the inductance per unit length along the current direction of the 
top plate. Since the electrodes have different widths and penetration depths, in general, 
Ct 7^ Cb- In Section III we will show that for our high quality all- Niobium LJTJs, having 
the base electrode thinner and wider than the top one, we found Ct ~ 3£fc. According 
to the theory of the two-conductor transmission lines^^, the inductance per unit length, 
£j, of a LJTJ, seen as a transmission line structure, is simply obtained as the sum of the 
inductances/unit lengths of the bottom and top stripes, i.e., 

/:.^A + A^.„(^ + |). (4) 

Historically, the boundary conditions for Eq.(|3]) were derived under the implicit assumption 
that Wb = Wt = W, so thalpIES] £^ _ jj^^^d^/w. However these conditions are not fulfilled in 
real samples, especially for window-type LJTJs used nowadays whose electrodes have quite 
different widths; typically, Wb > Wt > W. 

The paper is organized in the following way. In Section II we will overcome this limitation 
by extending the existing theoretical modeP° to LJTJs having different electrodes widths, 
strictly speaking, different inductances per unit length. At the same time we will derive 
the most general boundary conditions for Eq.(|3]) needed to correctly describe any self-field 
effect in LJTJs. Next we will focus on the specific case of LJTJs with doubly connected 
electrode(s). Later on we will consider the consequences of the fiuxoid quantization and 
energy minimization principles. In the next section we will describe our experimental setup 
and our samples; in addition we will present their magnetic diffraction patterns and discuss 
how the experimental findings can be unambiguously interpreted in term of our modeling. 
Finally, the conclusions will be drawn in Section IV. 

II. THEORY 

Figur^a) displays the top view of an in-line LJTJ having the most general biasing config- 
uration. Only three of the four currents Jj indicated in the figure are independent, since the 
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(a) (b) 

FIG. 2. (a) Top view (not to scale) of an in-line long Josephson tunnel junction in the most general 
bias configuration, (b) Cross section (not to scale) of the junction area showing the currents in the 
electrodes. The base electrode is in black, the top electrode in gray and the tunneling insulating 
layer is white in (a) while it has a wavy pattern in (b). 

charge conservation requires that Ii + 1^^ = I2 + I^- In addition, the /j's can be expressed in 
terms of the net current / crossing the tunnel barrier as I = Ii — 12 = 1^ — 14^. The junction 
cross section is sketched in Fig|2](b) together with the current distribution at the input and 
output end of the Josephson structure and along the junction electrodes. Here, Ib{,X) and 
It{X) denote, respectively, the local supercurrent flowing parallel to the insulating layer 
in the bottom and top junction electrodes with X G [— L/2,L/2], so that Ji = /(,(— L/2), 
I2 = Ib^L/^), I3 = /f(L/2), and = /<(— L/2). Next we observe that, due to the charge 
conservation, the total current, /b(X) + It{X), through any junction cross section in the 
Y — Z plane must be constant. For in-line LJTJs, it is important to distinguish between 
the symmetric and fully asymmetric configurations: in the former, the bias current / enters 
at one extremity and exits at the othei'^^'^J^ (/ = ±Ji = iJa and I2 = h = 0), while in 
the latter, the bias current enters and exits from the same 

extremit}^™!!] = = ^/^ 

and I2 = I3 = 0) . We will analyze the general cases when all /j 7^ 0. The coordinate system 
used in this work is indicated in Figs. |2](a) and[2]^b). 

A. Boundary conditions 

Figure |3] shows one elementary cell of the equivalent circuit for a static Josephson junc- 
tion transmission line. Classically, the two inductances have been merged in their parallel 
combinatioEp2Ell so that the role played by each supercurrent separately was lost; however, 
for our purposes it is mandatory to keep the distinction, since, in general, L^i^ C-t- In the 
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absence of an external in-plane field, Happ = 0, the magnetic flux A$j linked to the cell is: 



A$j = CtAXIt - CbAXIi 
wher^ $j = $00/271. Then, in the limit AX 0, 

By differentiating Eq. in force of Eq. (|3]), we end up with: 



Interestingly, by integrating back the equation above over the junction length L and consid- 
ering that, according to the notations of Fig. |2| Jf(L/2) = J3, Ji(— L/2) = I^, Ih(L/2) = I2, 
and /(,(-L/2) = Ji, we find Cb+Ct = fiode/W, i.e. \b/Wb+\t/Wt ~ {\b + \t)/W; clearly, this 
can only be acceptable if W ^ Wb ^ Wt. To overcome this limitation, we introduce a new 
(smaller) effective barrier penetration d'^ = WCj/hq that takes into account the screening 
effect expected when the junction electrodes are wider than the tunneling barrier, namely: 

de = Afc— + At— = \ , (6) 

Wb Wt Wb wt 

where Wb,t denotes the width ratio Wb^t/W. Being Wb^t > 1, then d'^ < \b + \t = de- We stress 
that a smaller magnetic penetration results in a smaller magnetic flux through the Josephson 
barrier (in small junctions or at the extremities of LJTJs), but not to a reduced amplitude 
of the magnetic field threading the barrier. Indeed, we believe that, due to demagnetization 
effects in the finite-thickness films, the amplitude of the magnetic field threading the barrier 
is larger than that of the applied field. The effects of a reduced magnetic thickness and 
an increased field amplitude partially compensate each other; however, this is not a good 
reason to ignore them. Since the screening and demagnetization effects depends, respectively, 
on the film widths and thicknesses, they are independent; we leave the investigation of 
demagnetization in Josephson structures to a future study. In the rest of this Section 
we will carry out our analysis substituting de with d'^ in the magnetic Josephson equation 
Eq. ([2]) ; consequently, this new magnetic thickness also enters in expression for the Josephson 
penetration depth and corresponds to the Aj inflation, already investigated using different 
approache^^^l*^, occurring in window-type Josephson tunnel junctions. For future purposes, 
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FIG. 3. Elementary cell of the equivalent lumped circuit for the static Josephson transmission line. 

we also introduce the two relative inductances per unit lengths A;, = Ch/Cj and At = 
Ctl Cj = 1 — and we note that, when Ch = Ct, then A;, = A^ = 1/2 as it was implicitly 
assumed in all previous analytical works on LJTJs. We anticipate here that for our samples 
we found quite different relative inductance per unit lengths, namely A;, ^ 0.26 and A^ ^ 
0.74. A practical expression for computing the bottom relative inductance not involving the 
junction width W is: 



XbWt' 



(7) 



Inserting Eq.([5]) into the magnetic Josephson equation ([2]), we obtain the local magnetic 
field in the barrier plane Hy{X) in terms of lb and If. 



Hy{X) = K- 



AtltiX) ~ AbhiX) 



(8) 



dX W 

Even when the junction electrodes are made of the same material and have the same thickness 
and quality, meaning that At = Xb, then the dependence on the electrode widths remains: 

^ hiX)Wb - IbiX)W, 

'~ w{Wb + Wt) • 

Equation ([s]) allows us to correctly derive the boundary conditions for the static sine-Gordon 
equation in Eq. ([s]) in the general case W ^ Wb and W ^ Wt and in the presence of an in- 
plane field Happ'. 



dX 



-L/2 



H, 



app 



Ath - Ab/i 
W 
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For Afe = A( = 1/2, we recover the boundary conditions by Owen and ScalapincP^ (/i = Is = 
/ and I2 = I4 = 0) that were generally adopted thereafter for untrue symmetry reasons. 
From the boundary conditions above, it follows that: 

where the last term, vanishing when = At and J4 + J3 = Ji + J2, has been omitted in all 
previous analysis of in-line LJTJs. We like to point out that Eqs.([9]) are very general and 
should be used to correctly describe the so-called self-field effects occurring in LJTJ. They 
also apply to LJTJs with mixed in-line and overlap biasing. Unfortunately, their imple- 
mentation requires the separate knowledge of the bottom and top electrode inductances per 
unit length (rather than their sum). Indeed, the inductance per unit length was analytically 
derived by Chan^^QlfQ]; ^ superconducting strip transmission line, i.e., a structure consisting 
of a finite-width superconducting strip over an infinite superconducting ground plane, as far 
as the strip linewidth Wt exceeds about the insulation thickness tox- Definitely his results 
can be used when Wb » Wt, but, unfortunately, no analytical expression is available when 
both electrodes have finite and comparable widths. 



1. Single loop 

Let us choose that the currents are positive when they flow from the left to the right junc- 
tion ends, i.e., counterclockwise in the case of a ring-shaped electrode sketched in Fig. [l|a). 
Then the boundary conditions for the single loop case can be derived as follows. With 
reference to Figs|2| we recognize that for the single loop configuration Ii = —Idr + o:I, 
I2 = —Idr ~ (1 ~ h = I and J4 = 0. Then Eqs.(|9| become: 



K 



dX 



Hy{-L/2) 



X=-L/2 



rr , ^(T _ n - Pf _ T- 

^app ' Oil) -Og ) 



K 



dX 



Hy{\./2) = Hy{-L/2) 



X=L/2 



w 



aAh 



(10) 
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with He = Happ + Hrad and 



Hrad = (11) 

Of course, if the loop is formed by the top, rather than the bottom, electrode. At should 
replace A;, in the above expression. Next, in this specific case, in order to have /;,(— L/2) = 
-Icir + al, /6(L/2) = -Icir - (1 - a)I, /t(-L/2) = 0, and /t(L/2) = /, it must be: 



hiX) = + al- W, [ .UX')dX' (12) 

J -L/2 

and 



It{X) = Wj I .UX')dX'. (13) 

J -L/2 

In normalized units of x = X/\j, the differential equation Eq.([3]) becomes: 

^=sin0(x), (14) 
with X G [—i/2,i/2] and i = L/Xj is the junction normalized length. Further, we will 



normalize the magnetic fields to JcAj, so that the boundary conditions Eqs.(lO) for a LJTJ 
with a doubly connected base electrode are: 



dx 



= hi = he - aAbi; — 

x=-l/2 '^■^ 



= hr = hi + i, (15) 

x=e/2 



where the term l = I / Iq = hr — hi is the external bias current / normalized to Jq = JcWXj. 
With these notations, the normalized critical magnetic field he of a short Josephson junction 
is 27i/i; further, defining i^ir = Idr/h and Kad = Hrad/Jc>^j = Abidr, then he = happ+Kad = 
happ + Ahieir- Wc uotc that what matters now is the product aA;,, rather than a itself and 
that the symmetry condition now corresponds to 2aAb = 1, which can never be achieved if 
Afe < 1/2. In the early eightie^, the reported asymmetric behavior of samples that were 
believed to be symmetric led many experimentalists to abandon the in-line geometry in favor 
of the overlap one. 
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2. Double loop 



We now consider the most general case of a LJTJ having both electrodes doubly con- 
nected. For the sake of simplicity, we now assume the two loops to be rectangular, as 
depicted in Fig. [l](b) where /3 is the fraction of the current / diverted in the left arm of 
the top-electrode loop (obviously, it is impossible to realize a topologically equivalent lay- 
out with two rings). As indicated, the magnetic fluxes $e.t and $e,fe linked, respectively, 
to the top and bottom loops induce the clockwise circulating currents Idr^t and Icir,b in the 
respective loops. Further, we recognize that Ji = —Icir,b + ctl, h = —Icir,h ~ (1 ~ «)-^) 
J3 = —Icir,t + (1 ~ (3)( and J4 = —Icir,t — (^I- From Eqs.([9]) the boundary conditions are: 



d4) 



= Hy{-L/2) 

X=~L/2 



a app + -^-^{Icir,t - (31) + -^-^{Icir,b " Oil) 



(3 At + aAb ^ 

lip j; 

W 



= Hy{h/2)=Hy{-h/2) + ^ (16) 

X=L/2 

with He = Happ + Hrad and Hrad = {Ablcir,b + Atlcir,t) /W: we note that the two circulat- 
ing currents interfere constructively since they flow in opposite directions, but also on the 
opposite sides of the tunnel barrier. The single loop configuration can be considered as a 
particular case of the double loop configuration in which /3 = (or 1). Similarly the free 



junction can be recovered by setting both a and (3 in Eq.(16) to any of their extreme values. 
In the rest of the paper we will limit our analysis to devices with the single loop configuration 
for which experimental data are available. 



B. Approximate solutions for LJTJs 

With the assumption that the Josephson junction is so long that the magnetic field in 



its center can be neglected, (j)x{0) ~ 0, an approximate solution of Eq.(14) is given by 
the superposition of two static non-interacting fractional fluxons pinned at the junctions 
extremitied^ 
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0(x) = 0'(x) + 0'-(x), (17) 

with (f)\x) = 4 sgnhi ta.ia'^ exp — (a; + + i/2) and (f)'^{x) = 4 sgn/i^ tan^"*^ exp(x — — ^/2), 
and sgn is the signum function. We observe that 0' and 0'' do not overlap, 0'(x) 0''(x) ~ 0, 
so that 0^(x) ~ [0'(x)]^ + [0^(a;)]^ (similar arguments hold for the phase derivatives). As an 
example, for i = 10, both |0'(x) 0''(x)| and |0^(a;) 0^(a;)| are everywhere less than 8 x 10""^. 
The phase profile in Eq.(17) can also be cast in the formP: 



sin = sgnhi sech (^ + ^' + 2] ~^ ^S^^r sech ( x + ^ 



From the phase derivative: 



(f'xix) = —2 sgn/i; sech ^ ^ 2 j ~'~ ^ sgn/i^ sech ^ — ^ 

we infer that C,i and are two non-negative independent constants set by the magnetic field 
at the boundaries hj./. 

hr,i = 4>x (=^0 = ^(^^/'^) ^ 2 sgnhr,i sech^^^z, (18) 

i.e., ^r,i = cosh~^ |2//ir ;|. This indicates that, in the Meissner state, the largest possible 
amplitudes of the boundary fields are- - hj. = hi = 2, corresponding to i = 0; then the 
junction critical field is he = 2, corresponding to He = 2JcXj. 



It is worth to remark that the solution in Eq.(17) only depends on the specific boundary 
conditions imposed by the system geometry. However, in the limit of vanishingly small bias 
current, the self-field effects disappear and the junction geometrical configuration does not 
affect the phase profile; in other words, for h^. = ±2, in-line, overlap and 5-biasecP^ LJTJs 



all have the phase profile given by Eq.(17) with hr = hi = ±2 



It can be easily proved that 0(0) ~ 0^(0) — and that, with ^i^r = 0, then 0(±£/2) 



sgnhr^iTT, meaning that Eq.(17) corresponds to a semifiuxon (vr jump) at each junction 
end, as shown in Fig. |4]^a) for i = 10 and hr = —hi = 2. For the second derivative 
of the phase we have: (pxx{x) = sin0 = — 2 sgn/i; sech(x + ^/ + £/2)tanh(x + ^« + 
+2 sgnhi sech(x - - ^/2) tanh(x + - ^/2). 



As \hi\ {\hr\) exceeds h^ then ^i {^r) becomes negative, the solution in Eq.(17) is no longer 



stable and we exit the Meissner regime; in fact, for C,i^r < the phase at the extremities grows 
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(a) (b) (c) (d) 

FIG. 4. (a) Phase profile (/>(x) of a long = 10) in-line symmetric {a = 1) Josepfison tunnel 

—hi = 2 that is = 0; (b) its first derivative 4>x, (c) its second 



junction as in Eq.(17) for hr - 
derivative (pxx, and (d) its normalized energy density TisC- 

above the threshold value = 7c/2 and in a dynamic scenario one (or more) integer vortex 
(fiuxon, antifiuxon) gradually develops at each extremity and moves toward the center under 
the effect of the Lorentz force so that some magnetic flux enters into the junction interior. 
Furthermore, for generic non-negative ^r,i values, the phase difference across the junction 
length A0 = 0(£/2) - (j){-£/2) is: 



2 Sin- I 



2 sin 



.1 hd 
2' 



(19) 



and corresponds to what has been called a /c-fractional vortex in RefJ ^^ * ^^ l, where k = A(/)/27r. 
Presently, semi {C,r.i = 0) and fractional {C,r,i > 0) vortices are receiving a great deal of 
attention in the context of — vr transition Josephson junctiond^^^. 



C. Junction energy 

By applying a Lagrangian formalismP^ to the sine-Gordon equation, it is possible to 
derive that the (static) Hamiltonian density HsG of a LJTJ is: 

^.G(a;) = y + l-cos<^, (20) 

in which the first term accounts for the magnetic energy stored in and between the junction 
electrodes and (1 — cos0) is the Josephson energy density associated with the Cooper-pair 
tunneling current. 



The phase solution, Eq.(17), for LJTJs also satisfies the equality: 2(1 — cos 0) = (pl, meaning 



that for very long junctions the Josephson energy density equals the magnetic energy density 
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FIG. 5. Numerically computed zero-field normalized junction energy HsG vs. of a symmetric 
(2aA{, = 1) in-line Josephson tunnel junctions with for different normalized lengths L For £ — t- oo 



the analytical expression in Eq.(23) was plotted. l = I/Iq is the normalized bias current. 



(this is not true for intermediate length junctions with i < 27r). Inserting the expression of 



the phase profile Eq.(17) in Eq.(20), then the junction energy density reduces to: 



-Hscix) = [Mx)]^ = 4sech2 { x + + 



4 sech^ \ X + 



(21) 



which is shown in Fig. Hid) for i = 10 and 



2 J V 2^ 

0. The junction energy HsG can be computed 



from Eq.(21 ) as a function of the boundary conditions h^r < 0, i.e., of the external magnetic 



field he and bias current t,: 



(22) 

where tanh^^,, = ^1 — sec\^{E,i.r) = a/1 — (/i«,r/2)^ and Eq = is the well-known fluxon 
rest mas^ normalized to Eq = /o$o/2vr that, depending on the junction's electrical and 
geometrical parameters, represents its characteristic energy unit; typically Eq is in the 10~^^J 
range, that is, several orders of magnitude larger than the thermal energy at cryogenic 
temperatures. In real units, the junction energy is Ej = EqHsg- 

For the particular case of a symmetric junction (2q;A{, = 1) in zero field {he = 0), then 



hr = —hi = l/2 so that Eq.(22) reduces to: 
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= 4(^/4)^ + {i/Af + 0.5(^4)' + 0{n. (23) 

Figure [s] displays the numerically computed zero-field energy Hgc versus for symmetric 
in-line junctions having different normalized lengths for £ > 27r the numerical findings 



are very well approximated by Eq.(23). As the expansion in the r.h.s. of Eq.(23) indicates, 
a LJTJ can be thought of as a non-linear inductance 2Ej/P which, in contrast to the 
small junction case, does not diverge at the critical current (this is because Ic < JcWh). 
The largest inductance value, achieved when the junction is biased at the critical current 
h = h,max = 4/o, is Lo = 2H sg{I c,max) / llmax = ^o/^o = ^of^nlo = CjXj and is the 
junction natural unit (typically a fraction of a picohenry) that will be used later on to 
normalize inductances. Normalizing the magnetic fluxes to the magnetic flux quantum $0) 
then the normalized circulating current can be written as idr = 27c{n — 4'e)/hoop, where 
(f)e = and lioop = Lioop/ Lq. With such notations, lb = Lb/Lo = Abi. 

Neglecting the mutual inductance effects, the system total energy Etot consists of two inde- 
pendent contributions Efot = Em + Hgc', the former is the magneto-static energy, Em, stored 
in the inductances Li and L2: 2Em = Lilf + L2/I; the latter is the previously discussed 
junction energy, Hsg, which, as said before, takes into account both Josephson energies and 
the magnetic energy associated with the bias current / flowing in the junction electrodes. 
In terms of normalized quantities, considering that Ii = al — Idr and I2 = —(1 — a) / — Idr, 
we have that Era = E^/Eq is: 

Em = [aHi + (1 - a)%] - + 



+ [(1 - a)k - ah] Li,„ + {h + h) '-f. (24) 
Em is minimum for a = amin, where: 



h '^cir h ^cir 

amin ~j I J r z z I , 

<1 + <2 I. lioop ~ l-b I' 

with L 0. In passing we observe that, in the absence of circulating currents, amin is 
independent of t. Moreover, in many cases of practical interest, Em >> 8 > Hsg so that 
the junction energy can be neglected. 
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FIG. 6. Magnetic diffraction patterns ic{he) for a very long in-line Josephson junction with different 
values of the symmetry parameter < a < 1 and for three Af, values, namely, 0.25, 0.5 and 0.75. 
The critical current is normalized to JcWXj and the magnetic field to JcAj. 

D. Magnetic diffraction patterns 



Setting hi and hr at their extreme values ±2 in Eqs.(15), we obtain the magnetic diffrac- 
tion pattern (MDP) ic{he) in the Meissner regime as a function of the symmetry parameter 
a: 



2 + he 

icihe) = — ^ for -2 <he< hmax, 

2- he 

for hmax <he<2, (25) 



1 — aAh 



with hmax = 2{2aAb — 1) being the field value which, for a given a, yields the maximum 
critical current idhmax) = 4. This v value can be achieved only when, as depicted in 
Figs. |4| /ir = —hi = 2. The h^ax expression turns out to be very useful in the experiments 
to determine the product aAf, from the analysis of the junction MDP: 

2 He 

in which the plus sign has to be chosen when h^ax is negative {2aAi, < 1) and vice versa. 
Figs. |6](a-c) display the MDPs for different values of the symmetry parameter < a < 1 
and for three Af, values, namely, 0.25, 0.5 and 0.75. As already stated, in general, he results 
from the sum of two contributions, namely happ and hrad- However, here we are neglecting 
the contribution to 0e (and hence to hrad) deriving from the current It running in the top 
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junction electrode; in the experiments, when needed, this contribution can be compensated 
by an external field perpendicular to the loop plane. 

We remark that the Zc(^e) patterns are piecewise linear and, in general, they have two quite 
different absolute slopes \dic/dhe\ on the left and right branches. The wide range of linearity 
is very attractive for the realization of cryogenic magnetic sensors with a large dynamic range 
especially because large slopes can be achieved in the [—2, hmax] interval. If the in-plane 
modulating field He is the radial field Hrad induced by a persistent current Idr circulating 
in the bottom doubly connected electrode, then H^, = Hrad = ^bldr/W. It is then possible 
to define a current gaiiP^^ gi = \dlc/dlcir\ = \dic/dicir\ = Af,\dic/dhrad\ = Ab\dic/dhe\, with 
icu- = Icir/lQ and hrad = Hrad/ Jc^j = ^bicir] that is: 



- for - 2 < /le < hmax 

a 



for hm.ax <he<2. (27) 



1 - aA 



In zero external field, varying a in the range [0,1], we have current gains in the range 
[Af,, Ab/ At]. Since in all real samples Wt < Wb, in the case of electrodes having the same 
effective penetration depth A, then At > Ab, suggesting that, for a given a, it is preferable 
to realize the loop with the top electrode, rather than with the bottom one; this can also 



be inferred by comparing Figs. [6] (a) and (c). Further, according to Eq.(27), large current 
gains can be achieved with small-a samples by flux biasing the loop to have h^. < h^ax- 

For each trapped fluxoid the currents circulating around the loop change by an amount 
Alcir = ^o/Lioop corresponding to a jump in the critical current: 



AJ, = g^AI^r = (28) 

^loop 

Each trapped flux quantum results in a small but detectable change in the junction critical 
current Ic- It becomes therefore possible to readout the number n of flux quanta trapped in 
a superconducting loop by means of an in-line LJTJ. It is worth to stress that the presented 
findings constitute an improvement in the state-of-the-art of current or magnetic sensorP^l. 
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E. Fluxoid quantization 



The internal magnetic flux $j within the loop is the sum of externally applied flux $e 
and the self-flux, $s = Lioophir, produced by the shielding current, Idr, which circulates 
around the loop to restore the initial flux: 

Also $e results from the sum of two contributions: 

$e = $a + $t = floHzAeff + Ct / It{X)dX. (29) 

J-L/2 

The first term is the applied or geometrical flux, due to the uniform magnetic field Hz 
externally applied in the direction perpendicular to the loop plane and ^e// is the effective 
flux capture area of the loop. For narrow loops, the pick-up areas can be well approximated 
by their inner areas. The second term, $t, is the non-linear flux contribution due to the 
non-uniform current It flowing in the junction top electrode. Further, the reaction flux, <l>s, 
can be also expressed in terms of the currents in three inductive paths of the loop, Li, L2 
and Lh'. 



= Lih + L2/2 + 



L/2 

(Li + L2)(a/ - Icir) - L2I + Cb I h{X)dX 

-L/2 



(30) 



where $6 is the non-linear flux contribution due to the currents It flowing in the junction 
bottom electrodes. 

In the presence of an external magnetic field, Happ, applied in the loop plane and perpen- 



dicular to the long junction dimension, L, by using Eqs.(12) and (13), after some simple 
algebra we end up with: 



(31) 



5* 

At 



CbL {al - Ic 



A(f) 
zn 



(32) 
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L /2 

with, as before, A0 = j_^^^(j)x{X)dX = 0(L/2) — 0(— L/2) being the Josephson phase 
difference across the junction and $' = ^odeLHapp a factitious flux threading the LJTJ 
barrier. 

The single-valuedness of the phase of the superconducting wave function around the loop 
(fluxoid quantization) required^ 



= + (Li + L2) (a J - hir) - L2I + $ J = n<l>o, (33) 

in which n is an integer number, called the winding number, corresponding to the number 
of flux quanta trapped in the ring and <l>j = + $5 is the non-linear contribution to the 
internal flux due to the currents It and lb flowing in the junction electrodes: 



/•L/2 /•L/2 

^j = Ct It{X)dX + Cb / h{X)dX 

J-L/2 J-L/2 



2AiA,L£j(a/ - /,,,) + (A^ - A;,) ( |^A0 - $' ) . (34) 



Inserting the expression Eq.(34) for $j into Eq.(33) and switching to normalized units, the 



fluxoid quantization law reads: 



AAA0 = {I2 - alioop)i + lbAA{aL - i^rr) + 27rAA0', (35) 

where AA = A^ — A?,. By its deflnition, |AA| < 1. If the material, quality and thickness of 
the junction base and top fllms are similar, then Xt ^ Af, and, consequently, the expression 
for AA reduces to AA ~ (Wb — Wt) / {Wb + Wt) > 0. If At = Aj, then, for symmetry reasons, 
the non-linearity of the Josephson element does not play any role, as far as concerns the 
fluxoid quantization and we flnd that, as expected, a is given by the ratio: 

— -^2 k /„„x 
"0 = J = -J — • (36) 

-'-'loop Hoop 

It is interesting to observe that the same result is obtained by minimizing the magnetic 



energy Em in Eq.(24), if li were replaced by h + lb = hoop — h- In other words, for AA = 0, 



both the fluxoid quantization and the energy minimization carry the same information, 



although they are independent principles. Eq.(35) is the reason why in our analysis we need 



to distinguish Cb from Ct so that, in general, the fluxoid quantization imposes a constraint 
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on the Josephson phase difference A0 across the junction. Provided AA 7^ 0, Eq.(35) can 
be rearranged as: 



A0 = AA Hioopiao - a)L + Ad{aL - idr) + 2n^', (37) 

The fiuxoid quantization leaves the parameter a as a still unknown quantity that has to be 
determined by using the energy minimization principle. Since the system total energy Etot 
is the sum of two contributions, it might have more than one local minimum corresponding 
to states with different critical current Ic- We remark that the fiuxoid quantization and the 
energy minimization are independent principles and as such they must be satisfied simultane- 
ously. Summarizing, the analysis of a LJT J having a doubly connected electrode requires the 



self-consistent solution of the differential equation in Eq. ( 14 ) with the boundary conditions 



Eq.(15), the constraint Eq.(37) and the further requirement of energy minimization with re- 
spect to the parameter a that will depend on the system parameters l, happ, Kad, h, h, h, A^ 
and can be considered as a system degree of freedom. This complex task cannot be carried 
out analytically and, in general, one should use rather involved numerical methods. 



The bond imposed by Eq.(37) frustrates the Josephson phase 0(x) along the junction. As 
far as l << 4, (f){x) can adapt its profile to the constraint. However, increasing l one 
reaches a point when no phase profile is compatible with the corresponding bond and a 
premature switching occurs. The capacity of handling the phase frustration grows with 
the normalized length of the junction. Ultimately, in the system under investigation, the 
fiuxoid quantization results in a phase frustration which, in turn, reduces the junction critical 
current. With a fixed phase difference, dynamic processes such as the resonant fiuxon 
motiori^ cannot exist any longer; nevertheless, flux flow processes will survive and their 
stability could even be enhanced by a fixed A0. 

For the double loop devices a and /3 can be considered as two degrees of freedom for the 
system and we shall have to apply two fiuxoid quantization rules which act as a double 
constraint on the phase difference A0 at the junction extremities. This is equivalent to fix 
the phases 0(— £/2) and 0(^/2) at the junction extremities. As shown in Fig. |4]^c), the peaks 
of the supercurrent density are not at the edges of the junction but are positioned where is 
needed to maximize the total current. Therefore a LJT J can carry a net supercurrent despite 
the boundary phase coercions; however, in this case, only internal dynamic processes, such 
as the resonant plasma oscillations, will be allowed. The two quantization rules together 
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with the energy minimization condition allow to determine both a and /3. However, with 
two loops, the system total energy should also include the mutual magnetic interaction. We 
postpone the thorough analysis of such devices to a future work. 



F. Remarks 



In our model of window-type LJTJs we have neglected the thickness of the oxide layer 
which is correct as far as we deal with the tunneling region. A more realistic picture of real 
devices should consider that in the idle region surrounding the tunnel area the insulation 
between the bottom and top electrode is provided by an oxide layer typically made of a 
deposited SiO^ layer and/or a anodic oxide. The total thickness of this layer is comparable 
or even larger than the electrode penetration depths A;, ^ (and might also be comparable 
with the strip width). In Ref.l^, each electrode was modeled as a parallel combination of 
two stripes having quite different oxide thicknesses resulting in a rather involved expression 
of the effective magnetic thickness. The expression of d'^ proposed in Eq. ^ should therefore 
be considered as just a first approximation which needs to be refined. Both magneto-static 
simulations and/or properly devised experiments could improve our knowledge on this topic. 
In this Section the consequences of the fluxoid quantization were derived for a LJTJ built on 



a narrow superconducting loop and the resulting phase constraint in Eq.(37) also contains 
a term proportional to the junction normalized length i. However, it is not clear which 
junction length maximizes (or minimizes) the effects of phase frustration. Further, our 
results cannot be extrapolated to the limit of small junctions which are not affected by 
self-field effects, i.e., hr = hi. In previous (not pubhshed) experiments, no evidence of the 
fluxoid quantization was observed in the magnetic diffraction patterns of small Josephson 
junctions built on large inductance loops h « hoop- Nevertheless, we expect a different 
behavior in the limit ^ hoop-, that is, when the junction energy dominates the loop energy. 



III. EXPERIMENTS 

A. Experimental setup 

Our setup consisted of a cryoprobe inserted vertically in a commercial LHe dewar (T = 
4.2K). The cryoprobe was magnetically shielded by means of two concentric Ph cans and 
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a cryoperm one; in addition, the measurements were carried out in an rf-shielded room. 
The external magnetic field could be applied both in the chip plane or in the orthogonal 
direction. In fact, the chip was positioned in the center of a long superconducting cylin- 
drical solenoid whose axis was along the y-direction [see Fig. |2]^a)] to provide an in-plane 
magnetic field Happ = Hy- A transverse magnetic field, Hz, was applied by means of a 
superconducting cylindrical coil with its axis oriented along the Z-direction; this trans- 
verse field induces a controllable shielding current Idr circulating in the superconducting 
loop that, in turn, generates a radial field, Hrad, in the insulating layer of the Josephson 
structure that algebraically adds to Happ. The field-to-current ratio was 3.9 /iT/mA for the 
solenoid and 4.4yuT/mA for the coil. These values have been numerically obtained from 
Comsol Multiphysics^^ magneto-static simulations in order to take into account the strong 
correction to the free-space solution due to the presence of the close fitting superconducting 
shielcpSl. The effects of a transverse field on the static properties of both shortP^ and lon^^ 
Josephson tunnel junctions of various geometries have recently been investigated. 

B. Samples 

In Figs. [7](a) and (b) we report the two (topologically equivalent) geometrical configu- 
rations used for our experiments; the rectangular loop in Fig. [7]^a) has a mean perimeter 
approximately equal to 27r times the mean radius i? = 53 fim of the ring-shaped loop in 
Fig. ^h). Since the semi-angular length S = L/2R 1, to a first approximation, the curva- 
ture of the junction built on top of the ring can be ignored. In both cases each loop allocates 
two in-line junctions sharing the doubly connected base electrode; they can be biased sep- 
arately and under different bias configurations corresponding to different L2/L100P ratios, 
i.e., ao values. In addition, if the bias current is applied through their respective counter- 
electrodes (terminals 1 and 4 in Fig. [?]), the two junctions are series biased. Here, we will 
only present experimental data on single loop devices with just a single junction biased; 
the common biasing of both junctions has been intended for a different class of experiments 
aimed to improve the results of Ref.^ and will be the subject of future work. 
To correctly analyze our geometrical configurations, we have to introduce one more path 
inductance corresponding to the loop section acting as the base electrode for the passive 
LJTJ. Furthermore, if both junctions were biased simultaneously, the fiuxoid quantization 
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(a) (b) 

FIG. 7. Geometrical layouts (not to scale) of our samples: (a) rectangular and (b) circular loop. 
The semi- angular length of the junction is 5 ~ 1. 



would involve the phase differences across both junctions and the total system energy would 
include the two junction's energies, as well. However, if only one of the junctions is biased, 
then the model developed for the single loop can be fully adapted with the caveat that 
properly calculated contributions have to be added to the inductive paths Li and/or L2. 

We evaluated the loop inductances considering the loops as isolated superconducting narrow 
loops: the results are 90 and 240pH for the rectangular and circular loop, respectively. 
However, these are overestimated values, since the portions of the loop covered by the 
counter electrode have a lower inductance per unit length. The experimental data will 
provide more accurate values. 

In the experiments, we used high quality Nb/Al — Alox/Nb LJTJs fabricated on silicon 
substrates using the trilayer technique in which the Josephson junction is realized in a 
window opened in a 200 nm thick Si02 insulator layer. We measured a large number of 
junctions with width W = 1.5 fim and length L = 100 /im and Josephson current density 
Jc — 3.6kA/cm^, as measured in small junctions realized in the same fabrication batch. 
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The two junctions on a given loop only differ by the longitudinal idle regions which are 
50 fim for JJA and JJC and 1 fim for JJB and JJD. The nominal thicknesses and widths of 
the bottom and top electrodes were, respectively, df, = 100 nm, Wh = 6 /im, dt = 350 nm, 
and Wt = 4/xm , so that, assumin^^ XL,Nb = 90 nm, it is A;, = 45 nm, Xt = 85 nm and 
the effective magnetic penetration d'^, as given in Eq.(|6]), resulted to be d'^ ~ 43 nm. The 
same values in Eq.Q give A;, = 1 — At ~ 0.26 corresponding to Ct ~ 3Cb and AA ^ 0.48. 
The value of d'^, together with the critical current density given above, yields Aj ^ 12.7/im; 
henceforth, our junctions had a nominal normalized length L/Aj ^ 7.9 and we can treat them 
as long (L > 27rAj), one-dimensional (W < Aj) Josephson tunnel junctions. (By assuming 
that the magnetic penetration was = Ah + Aj ~ 130 nm, we would get A j ~ 7.3 /xm and 
Afe = \b/de = 1-At^ 0.35.) 

C. Magnetic diffraction patterns 

On real samples, the measurement of the maximum supercurrent versus the external field 
yields the envelop of the lobes, i.e., the current distribution switches automatically to the 




(a) (b) 

FIG. 8. Experimental magnetic diffraction patterns of a long Josephson tunnel junction built on 
top of a rectangular loop (sample JJB) under different bias configurations and applied magnetic 
field orientation: (a) in-plane and (b) transverse magnetic field. With reference to Fig. [7|(^b), the 
voltage is measured across the terminals 1 and 2, while the bias current is applied through the 
contacts: 1 and 2 for curves I, 1 and 3 for curves II and 1 and 5 for curves III. /q = 460 /iA, 
Hy,c = 490 A/m, and Hz,c = 80 A/m at T = 4.2 K. 
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mode which, for a given field, carries the largest supercurrent. Sometimes, for a given ap- 
plied field, multiple solutions are observed on a statistical basis by sweeping many times the 
junction current- volt age characteristic. Figs. |8](a)-(b) display, respectively, the in-plane and 
transverse MDPs for a LJTJ built on top of a rectangular loop under different bias con- 
figuration (resulting in different L1/L2 ratios). Qualitatively similar results (not reported) 
were obtained for samples with annular geometry. In both figures, the data labels I, II and 
III indicate that the junction bias current was applied through, respectively, the terminals 
1 and 2, 1 and 3 and 1 and 5. We remark that, as expected, all the MDPs are symmetric 
with respect to inversion of both the junction bias current, /c, and the applied magnetic 
fields, Hy or Hz- In Figs. |8[a)-(b) Jq = 460 /iA, Hy,c = 490 A/m, and Hz,c = 80 A/m. As 
shown by the dotted line, the critical field values were obtained by extrapolating to zero the 
linear branches of the MDP main lobe (by definition the critical field does not depend on the 
bias configuration). The theoretically expected value of Iq = JcWXj ~ 690 /xA exceeds its 
experimental counterpart by 50%. Since in the fully asymmetric configuration I the self-fiux 
effects are absent (with L2 = = 0, no fraction of the bias current circulates in Li), at a 
qualitative level, we recognize this discrepancy as the signature of the phase frustration in 



Eq.(37) provided by the fluxoid quantization. Put in a different way, the presumed harmless 
fact to have a doubly connected electrode depresses the critical current with about 35%. 
Unfortunately, even in absence of external fields and with the winding number n set to zero, 
no analytical expression can be found for the value of Iq expected in presence of the doubly 
connected electrode. 

Further, the expected value of the in-plane critical field, Happ,c = 2JcAj — 915 A/m is al- 
most twice as large as the measured Hy^c] this discrepancy can be ascribed to the previously 
mentioned in-plane demagnetization that, in the barrier proximity, squeezes the field lines of 
any magnetic field applied in the junction plane. Interestingly, good agreement was found 
instead for the transverse critical field, Hzci which makes us confident about demagnetiza- 



tion effects for the in-plane fields. In fact, according to Eq.(ll), the radial field experienced 
by a LJTJ built on a superconducting loop is: 

Tj _ ^bhir f^O^b^effHz , s 

W WLioop 

in which Lioopldr ~ /Uo-f^z^e// is the magnetic fiux threading the loop and Aeff the effective 



flux capture area of the loop. By deflnition also Hrad,c = '^Jc^j, then Eq.(38) provides the 
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following expression for the expected transverse critical field, if| 



^th ^ 2JcXjWLioop 
''^ /ioAfoAg// 

For narrow loops, as in our cases, the pick-up areas can be well approximated by their areas; 
for the rectangular loop of Fig. |8](a), ^e// ~ 4.5 x 10~^ mm^. By inserting the correct values 
in the last equation we get H^zc ~ 84A/m. According to Eq.(ll), when Hz = ^Hz,c then 



the shielding current is Idrc ~ ±2/o/A;,; at variance with the applied supercurrent /, the 
circulating currents are not bonded to the [— 4Jo,4/o] interval. 

We note that the experimental in-plane MDPs in Fig. [sj^a) closely reproduce the theoretical 
ones of Fig. [6|^a), apart from the fact that the largest critical current values Ic,max should 
be independent on the biasing terminals. However, the biasing configuration II and III are 
characterized by larger and larger circulating currents which depress the junction critical 
current even further. As already stated, our derivation of the MDPs did not take into 
account the self-flux effects. However, we were able to compensate these effects by a proper 
small transverse field Hz- 



According to Eq.(26), the product aoA;, can be determined from the field value H^ax corre- 
sponding to the largest critical current Ic,max- We stress that the position of H^ax remains 
independent on any possible suppression of the critical current Ic. Taking the fully asym- 
metric configuration I as a reference, ao = when H^ax = He, the obtained values of aoAj, 
for the biasing configurations II and HI are, respectively 0.11 and 0.21. With A;, ^ 0.26, 
we end up with ao = 0.42 and 0.81, respectively, for the biasing configurations II and HI, 
refiecting the fact that, as expected, the configuration II belongs to the asymmetric range, 
a < 0.5, while the configuration HI belongs to the symmetric one, a > 0.5. Furthermore, 
zero-field current singularities were observed in the junction current- voltage curves that re- 
mind of the resonant fiuxon motion; however, due to the phase torque induced by the fiuxoid 
quantization that mimic the effect of an external magnetic field, we believe that they are 
better ascribed to the asymmetric fiuxon propagation of Fiske-type resonances, observed in 
LJTJH3. 

We can now comment on the transverse MDPs in Fig. [8|b) in which drastic Ic changes 
are found but only for negative Hz values. In this case the externally induced circulating 
currents become increasingly important and, as discussed in the Section II, we have to 
consider the effects of the energy minimization with a being a degree of freedom bonded to 
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FIG. 9. Variation of the critical current as a function of the magnetic flux with which the loop was 
field-cooled through the superconducting transition temperature. Hy = Hz = and T = 4.2 K. 



the [0, 1] interval. It can be shown that, for generic hj- i values, the junction energy in Eq.(22) 
is minimum when 2aA = 1 + 2hrad/t' (we note that happ = Hy = in this situation). For 
concord hrad and l, the parameter a is squeezed to its upper value 1, while for hrad < —l/2 
it is bounded to zero; however, with the ratio 2hrad/i- G ["IjO]) the free parameter a will 
assume the value(s) in the range [0, 1] that minimize the energy. It is possible that, for a 
given hrad {Hz in the experiments), two values of l (J) exist that correspond to different 
junction energy levels. The coexisting states are evident in the transverse MDPs for for in 
the range [—0.5, —0.25]Hz^c- Therefore, we beheve that the sudden changes in Ic{Hz), that 
are absent in the in-plane MDPs, correspond to abrupt modifications of the phase profile 
4>{x) along the LJTJ and indicate that the energy minima correspond to quite different 
a values. Numerical simulations are planned to understand the experimental MDPs at a 
quantitative level. Finally, the appearance of displaced linear slopes for transverse magnetic 
fields larger than the critical value constitutes one more indication that the phase twist is 



increasing with the circulating currents, as suggested by Eq.(37) 



D. Signal-to- noise ratio 

With our devices a very large signal-to-noise ratio can be achieved in the detection of 
magnetic flux quanta trapped in the loop. Fig. |9] shows how the zero-field critical current 
changes of anyone of the two LJTJs on top of the rectangular loop when the system is cooled 
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through the NS transition in the presence of a transverse magnetic field which is incremented 
by steps corresponding to a small fraction of the magnetic flux quantum, A$e = 0.02 $o- 
Once the transverse fleld is removed, the quantized levels of the critical currents with A/^ ~ 
9.5 /iA are clearly visible out of the ±1 /xA rms current noise from the thermal fluctuation 
of the critical current. 

As the MDPs of Fig. [s] indicate, as far as \h(,\ < \hmax\, the slope of the positive and negative 
critical currents, respectively, and J~, is the same. Consequently, any circulating current 
Icir modulates and concordly, meaning that the offset current 1°^-^ = 1^ + 1^ changes 
twice faster. The root mean square noise of I"^^ is \/2 times larger than that of a single 
critical current, meaning that the signal-to- noise ratio is enhanced by a factor a/2. In 
addition to this, since the change with temperature of is numerically the same but 
opposite to that of J~, one more advantage of the current offset, as compared to just one 
critical current, is its much reduced sensitivity to any temperature drift that might occur 
during the measurements. The data shown in Fig. prefer to /°^-^/2. 



An accurate value of the loop inductance can be obtained from Eq.(28): 



-'loop 



AL 



For the rectangular loop with the bias conflguration I [gi = 0.33) we flnd Lioop = 72 pH. For 
the circular loop with the same biasing conflguration {g^ = 0.35), we measured A/^ = 5.2 /xA 
corresponding to Lioop = 140 pH. 

The superposition of the levels is due to the non-adiabaticity of the thermal transition^ 
and the transition from one level to the next follows a Gaussian probability law^®. In fact, 
during the normal-to-superconducting transition the loop temperature changed at a rate of 
5 X 10^ K/s. Indeed, this method is strongly inspired by the results found in our investigation 
of the spontaneous fluxoid formation in superconducting loops based on the detection of the 
persistent currents circulating around a hole in a superconducting fllm, when one or more 
fluxoids are trapped inside the hol^. 



IV. CONCLUSIONS 

In this paper we have revisited the theory of the self-fleld effects that characterize the long 
Josephson tunnel junctions and made them interesting for the investigation of non-linear 
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phenomena. Our analysis goes beyond the previous works in two ways: (i) it takes into 
account the different inductances per unit length of the electrodes forming the junction and 
(ii) it provides the boundary conditions for the most general junction biasing configuration. 
We applied the theory to the specific case of long Josephson tunnel junctions with not sim- 
ply connected electrodes. Apart from their intriguing physical properties, the interest for 
LJTJs built on a superconducting loop stems from the fact that they were successfully used 
to detect trapped fiuxoids in a cosmological experiment aimed to study the spontaneous de- 
fect production during the fast quenching of a superconducting loop through its normal-to 
superconducting transition temperatur^. We found that the single-valuedness of the phase 
01 of the order parameter of the bottom (or top) superconducting electrode (fluxoid quan- 
tization) gives raise to a variety of unexpected non-linear phenomena when coupled to the 
sine-Gordon equation for the phase difference 02 — (pi of the order parameters in the super- 
conductors on each side of the tunnel barrier. The principle of energy minimization was also 
invoked to determine the possible states of the system. We have focused on static phenom- 
ena such as the reduction of the junction critical current and its dependence on magnetic 
fields applied in and out of the loop plane. Nevertheless, also the dynamic properties, such 
as the propagation of non-linear waves, are expected to be drastically affected by the doubly 
connected electrode(s). Our experiments unambiguously corroborate the analytical findings 
and provide hints to implement the modeling. Future work should go in the direction of 
investigating the consequences of the fluxoid quantization for small and intermediate length 
Josephson tunnel junctions and on the dynamic properties of long junctions. 
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